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In response to a recent manuscript [lj on the analysis of interference patterns produced by ul- 
tracold atoms released from an optical lattice, we point out that in the presence of a Bose-Einstein 
condensate the interference pattern can be strongly modified by interaction effects and the presence 
of a harmonic trap superimposed on the lattice potential. Our results show that the visibility of 
the interference pattern is significant only if a sizeable condensate fraction is present in the trap, in 
strong contrast to the findings of ref. [J. 
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Introduction- Experiments with ultracold quantum 
gases in optical lattices rely heavily on time-of-flight 
expansion to probe the spatial coherence properties of 
the trapped gas @, 0, 0, 0] . When the phase coherence 
length is large compared to the lattice spacing, the post- 
expansion density distribution reveals a striking interfe- 
rence pattern with the same symmetry as the reciprocal 
lattice. As the phase coherence length decreases, e.g. on 
approaching the Mott insulator transition, the visibility 
of this interference pattern decreases accordingly, and for 
very short coherence length the underlying lattice order 
can only be revealed through higher-order correlations 
Q . A recent manuscript [H has raised the question, whe- 
ther even a minute condensate fraction on the order of 
l/(eN 2 / 3 ) (N being the total particle number) should 
give rise to an interference pattern with visibility V very 
close to unity, V ss 1 — e. According to the calculations 
in [![, a system of, e.g., N = 10 5 atoms should show an 
interference pattern with a visibility of V f=a> 99% for 
minute condensate fractions as 4% (ss 5000 atoms). As 
current experiments with ultracold quantum gases do not 
observe this (the best visibilities that have been achieved 
in our group being close to 90%), the authors conclude 
that the current experiments must be carried out in a 
thermal regime with no condensate or superfluid fraction 
present in the trap. 

In this paper we would like to point out that this re- 
verse conclusion is an artifact of the model used in [l[ to 
describe the time-of-flight expansion. A central assump- 
tion in Ref. [H is to consider that the condensate mo- 
mentum distribution is Heiscnberg-limited by the size of 
the system prior to expansion, and that this distribution 
is preserved in time-of-flight. Below, we would like to ex- 
plain that several effects spoil this assumption, as under 
realistic experimental conditions, one needs to take into 
account : 

- (i) The inhomogeneous trapping potential for the 
quantum gas in the lattice, 



- (ii) The mean-field broadening of the momentum 
peaks during time-of-flight expansion, 

- (in) The finite time-of-flight in the experiments. 
Practically, these effects tend to broaden the peaks 

present in the interference pattern, hence diminishing its 
peak value. In the presence of a significant non-condensed 
fraction, this will also lower the visibility. In cases when 
the initial momentum width is very large, for example 
for a cloud close to or in the Mott insulator regime, or a 
thermal gas well above the critical temperature, they can 
usually be neglected. However, the momentum width for 
a condensate is conversely very narrow. Even without an 
optical lattice, effects (i)-(iii) have been known to play 
an important role in the interpretation of time-of-flight 
images (see, e.g. 0,0, 0]). We point out that a number 
of other effects might be relevant in a given experimen- 
tal situation, for instance a finite imaging resolution, or 
collisions during the initial expansion phase that would 
scatter the atoms out of the condensate. Here we do not 
consider these aspects further, and consider an idealized 
experiment limited only by the three points above. Below 
we outline our arguments and point out the most critical 
steps to arrive at the conclusions given in We pro- 
pose a simple model to explain why the visibility of the 
interference pattern follows the condensate fraction in the 
system at least in a qualitative way. We also find that for 
the parameters that correspond to current experiments, 
the visibility of an expanding, non-condensed thermal 
cloud is always significantly lower than 1, a conclusion 
consistent with experiments 0, 0, 0, 0] . 

Setup and Optical Lattice Potentials - The optical 
lattice potential can be expressed as 

Voi» = V (sin 2 (/c L x) + sin 2 (fc L y) + sin 2 (/s L z)) . (1) 

Here Vq is the lattice depth expressed in units of the 
single-photon recoil energy E-r = h 2 /2mX 2 j . Here = 
27t/Al is the laser wavevector, Al is the laser wavelength 
and m is the atomic mass. In addition to the lattice po- 
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tential, an "external" potential is present, due to the pre- 
sence of a magnetic trap and also to the optical confi- 
nement provided by the Gaussian-shaped lattice beams. 
The external potential is nearly harmonic with trapping 
frequency 



A(2V - V%) 



(2) 



where u> m is the frequency of the magnetic trap, assumed 
isotropic, and where w is the waist (1/e 2 radius) of the 
lattice beams, assumed identical for all axes. This for- 
mula, slightly different from the one given in takes 
into account the modification of the vibrational ground 
state energy in each well due to the spatial variations 
of the laser intensities on the scale w (see 0), which 
amounts to a few percents change in wt- 

After releasing the cloud from the trap, the density 
distribution of the expanding cloud after a time of flight 
t is usually taken to be proportional to the momentum 
distribution, 



»(k)= \w (h)\ 2 s (k) , 



(3) 



i—x^y^z 



with a scaling factor r t . .f. = (fii/m)k. Eq. ([3]) is va- 
lid provided that two conditions are met : interactions 
have a negligible influence on the expansion, and the time 
of flight is long enough to treat the initial particle dis- 
tribution as point-like ("far-field" approximation). The 
envelope \w\ is the Fourier transform of the Wannier 
function in the lowest Bloch band, which we approxi- 
mate by a Gaussian, |u)(fc.j)| 2 w -^fr exp (— fcfiug), with 
wq the size of the on-site Wannier function. S(k) is gi- 
ven by <S(k) = J2i j e* k *^ ri_r;i ^(a|aj). Experimentally, one 
records the momentum distribution integrated over the 
probe direction z, which we call n±. We assume for sim- 
plicity that the width of S (k) is lower or comparable to 
1 /wq (this should hold in the situations where the popu- 
lation in the first excited Bloch band is negligible). Then, 
substituting \w (k z ) | with \w (0)| leads to 



(4) 



.(kj_) pa Ad / dk z S(k) 



where we introduced the lattice spacing d, a global factor 



a = ^- d \u,(k x )nw(k y ) 

[3| follows from 

n 

V = - 



The visibility as defined in 



(5) 



where n max = no(0) + nth(0) is measured at the first 
lateral peaks of the interference pattern, e.g. at K max = 
(1/2, 0) x 2ir/d, and where the minimum density n m i n = 
?ith(K m i n ) is measured along a diagonal with the same 
distance from the central peak, i.e. at the point defined 



the vector K r 



(1/V2,1/V2) x 2ir/d. 



Orders of magnitude - When condensed and non- 
condensed components coexist, the planar momen- 
tum distribution n± splits naturally into a condensate 
contribution n Q (r±) and a non-condensed contribution 
n nc (r±). Let us first consider the relative scaling of the 
condensate and non-condensate contributions. We write 
the condensate momentum distribution, characterized by 
the condensed atom number No and a momentum width 
Afco, as 



n (k ± ) = A 



No 



fo 



(Afc -d) 2JU \Ako) ' 



(6) 



where fo is a periodic, dimcnsionlcss function. Similarly, 
introducing the number of of non-condensed atoms N nc , 
we write the non-condensed contribution as 



71 n c(kj_) = A 



N, 



In 



(Ak nc -d) 2jnc \Ak 



(7) 



We assume the condensate contribution to produce a 
sharp interference pattern (/o(K max ) = 1 and /(K m ; n ) ps 
0), whereas the non-condensed part extends over the 
whole Brillouin zone (Afc nc w 1/d). The visibility then 
reads 



V = 



N + (AM) 2 A^nc / ( 



-) 



iVo + (Afc d) ^nc / ( + ) 

where we introduced 



(8) 



= (/nc(K max ) ± / nc (K mi „)) // (K max ) (9) 

which take values of order unity. The "suppression fac- 
tor" (Akod) 2 sets the importance of the non-condensed 
fraction ; if this ratio is very small, then even a tiny 
condensed number will yield a visibility close to unity. 

Uniform system - The system considered in [l| is a 
gas in a box with N s lattice sites in each direction (L 
in the notation of Ref. [H). The momentum width for 
the condensate is diffraction limited to Afco ~ l/N s d. 
Essentially this amounts to the assumption that the post- 
expansion momentum distribution of the condensate is a 
delta function. In [l| , any structure in the non-condensed 
component is neglected, so that / is uniform. Then, 



(Ak dY 



1 

TV 2 



10' 



(10) 



for a typical N s ~ 100. This is where the result in Eq. (4) 
of Ref. [l[ comes from, up to a factor wo/d which de- 
pends on how the k z integral is handled. In Fig. Q] (upper 
panel) the momentum distribution for such a case is dis- 
played. Essentially, the momentum peaks approach delta 
functions, explaining why tiny condensate fractions can 
lead to visibilities very close to unity. We note that a 
better way of analyzing the visibility, is to measure the 
weight, i.e. the integral over a peak over a finite area, 
which accounts for finite resolution effects that are inevi- 
tably present in experiments. 
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Fig. 1: Planar momentum distribution for fc H = 0. The upper 
panel is plotted for a uniform ideal gas with 7V S = 10 6 sites 
and one atom per site, and for a condensed fraction f c ~ 5 %. 
The lower panel is plotted for our harmonic trap model with 
interactions (using a = 0), also for a condensed fraction f c ~ 
5 % and a total number N = 10 6 using the trap parameters 
given in Note the vastly different scales on the vertical 
axes. 



Ideal gas in a harmonic trap - As mentioned above, 
current experiments are carried out in the presence of a 
harmonic trap superimposed on the lattice. In this case, 
one expects the ground state wavcfunction to be locali- 
zed on a scale given by the harmonic potential ground 
state extension ctt = \/h~/muJT, and not by the "num- 
ber of lattice sites" (roughly equivalent to the size uuq 
of the laser beams). In fact, usually one has wo 3> ctt- 
The main point here is that the suppression factor is no 
longer hopelessly small, especially when most atoms are 
non-condensed. However, as well known, this ideal gas 
model for the condensate is a poor description of realis- 
tic experimental situations 0|. 

Interacting gas in a harmonic trap - To obtain a bet- 
ter description, it is necessary to include interatomic in- 
teractions in the picture. In the Thomas-Fermi limit, 
which usually applies to experimental situations, the in- 
trap spatial condensate wavefunction will be significantly 
broader than in the ideal gas case, implying a narrower 
momentum distribution. However, as argued above, the 
post-expansion momentum profile can be very different 
from the sharply peaked in-trap profile. First, the inter- 



action energy is released as kinetic expansion energy, and 
second, for a finite time-of-flight t the initial size of the 
distribution is not necessarily negligible. Due to these ef- 
fects, the actual momentum distribution will differ from 
the one calculated using Eq. 

Here, we handle them in a heuristic way, by approxi- 
mating the maxima in the planar momentum distribution 
for the condensate as Gaussians, 



n Q (kj 



K„ 



N TT 2 A 
(Afceffd) 2 



exp(-ki/Afc c 2 ff ). (11) 



We assume an effective spatial width given by the quadra- 
tic sum of the initial width uo and of the expansion width 
HAko jmt corresponding to the release energy. In momen- 
tum space, this corresponds to Afc 2 ff = Afc 2 + ■ 
In the Thomas-Fermi limit, both the post-expansion mo- 
mentum width Afc and the initial size o§ arc determined 
by the interaction energy, which itself is proportional to 
the chemical potential fi. Accordingly, one has 



h 2 Ak 2 



fin, 



(12) 



where the coefficients a and fi are of order unity. The 
final result is 



Afc c 2 ff 



{oj T tf 



(13) 



corresponding to a suppression factor (Afc ef jd) propor- 

,2 

tional to Jr 2 ^ ~ h/Er < 1 much larger than expected 
from Eq. (|10[) (we estimate below the proportionality fac- 
tor). Hence, one expects the observed visibility to follow 
the evolution of the condensed fraction at least qualitati- 
vely. The argument also holds for the quantum-depleted 
part of the non-condensed cloud. 

Visibility for a trapped gas model To make the prece- 
ding arguments more precise, we introduce a semiclassi- 
cal description of the thermal cloud, for which interaction 
effects can usually be neglected to a first approximation 
The thermal component is characterized by a distri- 
bution function in phase space given by 



P(r,k) = 



cxp[(e k +V(r)-/i)/k B T]-l ' 



(14) 



with the tight-binding dispersion relation = 
2J^2 i=x y z (1 — cos(fcjd)), with the trap potential 
V(r) = ijTiw^r 2 . We can obtain the thermal fraction 
by integration over phase space, 



bh 



(2*)»/»(^) 8 £SLi 



nO-6J)/fc B T 



^I/ (2nJ/k B Tl\5) 



Here Iq is a Bessel function and i?th 



k B T 



the ther- 



mal cloud radius. The critical temperature for condensa- 
tion follows from Eq. fT5|) by setting the chemical poten- 
tial \i = Integrating the phase space distribution 
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over position and k z gives the planar momentum distri- 
bution, 



^ ' n—1 

I (2nJ/k B T)e~ n£k ± /kBT 



3 _oo_ e „( M _2J)/fc B T 



,3/2 



(16) 



This takes the form announced in Eq. ([8]). 

To estimate a and (3 in Eq. (|13p . we impose that our 
Gaussian model should reproduce the same interaction 
and potential energy as a condensate in the Thomas- 
Fermi limit, with chemical potential given by (Io| 

2/5 



-(. 
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Vl6\/27r 



AW 



Id 2 ) 



2/5 



(17) 



The potential energy per particle 0] i£, 



pot 



jmui^cr 2 gives a coefficient /3 ps 12/21. The coefficient 
a depends on how the interaction energy is redistribu- 
ted during the expansion. We consider the two limits, 
where the interaction energy is entirely concentrated ei- 
ther into the single diffraction peak under consideration 
(Eim = 2^/7 = 3h 2 Ak 2 /4m, or a = 8/21), or in the cen- 
tral peak avoided by the visibility measurement (a — 0). 
This provides upper and lower bounds for the actual vi- 
sibility. The results of the calculations is plotted in Fig. 
[2] for various assumptions : no trap, ideal gas in a trap, 
interacting gas in a trap with a = 0ora = 8/21. Com- 
paring the two latter cases, one sees that the final result 
depends sensitively on the value of a. A more complete 
calculation would probably yield a curve lying in between 
the two bounds plotted. Nevertheless, irrespective of the 
particular model, the curves demonstrate that the visibi- 
lity falls well below unity in a broad temperature range, 
where the condensate fraction stays non-zero. 

In conclusion, we have calculated the visibility of the 
interference pattern observed when releasing a finite tem- 
perature, Bose-condensed gas from an optical lattice. 
Using a realistic model for the condensate momentum 
distribution, we find that the evolution of the visibility 
with temperature or lattice depth follows the evolution 
of the condensed fraction at least qualitatively. We there- 
fore disagree with the conclusion made in Ref . [l[ , that a 
measured visibility deviating from unity necessarily im- 
plies the absence of a Bosc condensate in the system. 
This claim originates essentially from neglecting interac- 
tion effects on the gas expansion, which is a valid ap- 
proximation close or in the Mott insulator regime, but 
a poor one for a system including a sizeable condensate, 
with or without [9( an optical lattice. We do however 
agree with the authors of [l[ on the importance of finite 
temperature effects in current experiments, as implied 
by recent experimental [H, [II], [II] and theoretical work 
la . 14, 15, 3, 17 1, and certainly consider this topic to 



warrant further investigations. 

We would like to thank Martin Zwicrlein for useful 
comments on the manuscript. 
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Fig. 2: Upper plot : Calculated condensed fraction for an 
ideal Bose gas in a periodic potential, with and without addi- 
tional harmonic trap. Lower plot : visibility of the interference 
observed after release from the trap. The dash-dotted line cor- 
responds to an uniform ideal gas, as in the model of [l[, the 
solid line to a harmonically trapped ideal gas. The long and 
short-dashed lines correspond to two models of an interacting 
condensate, which should be seen as upper and lower bounds 
for the visibility (a — and a — 8/21, see text). The curves 
are plotted for N = 10 6 atoms, Vb = 10 Er and an expansion 
time t — 20 ms. 
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